Simulations of astronomical imaging phased arrays 
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Abstract 

We describe a theoreticai procedure for analyz- 
ing astronomical phased arrays with overlapping 
beams, and apply the procedure to simulate a 
simple example. We demonstrate the effect of 
overlapping beams on the number of degrees of 
freedom of the array, and on the ability of the ar- 
ray to recover a source. We show that the best 
images are obtained using overlapping beams, 
contrary to common practise, and show how the 
dynamic range of a phased array directly affects 
the image quality. 

1 Introduction 

Phased arrays have had many applications in 
the field of radar. In recent years there has 
been a considerable increase of interest in ra- 
dio astronomy applications, with such projects 
as the Square Kilometer Array (SKA) [HUT], the 
Electronic Multibeam Radio Astronomy Con- 
cept (EMBRACE) [2], the Karoo Array Tele- 
scope (KAT), and the Low Frequency Array 
(LOFAR) [9 . Furthermore, there are plans to 
place phased arrays in the focal or aperture 
planes of single dish telescopes, and as technol- 
ogy improves, phased arrays will be utilized at 
submillimeter wavelengths in addition to the ra- 
dio and microwave spectrum. 



In designing phased arrays it has been stan- 
dard practise to ensure that the synthesized 
beams are mathematically orthogonal, and thus 
when the incoming field matches one of the beam 
patterns precisely, only one port of the phased 
array is excited. Approximate (though not pre- 
cise) orthogonality of the beams is technically 
straightforward to achieve, and simplifies the 
analysis, but, as we shall show, is not necessarily 
desirable. In particular, in order to achieve or- 
thogonality one needs either to have the beams 
widely spaced, which leaves large parts of the 
sky uncovered, or to make the beams as narrow 
as possible, which increases the sidelobe levels. 
Recent work on shapelets [71 El IISI El E] and 
other non-orthogonal functions suggest on the 
contrary that it is desirable for the beams to 
overlap, such that is not possible for an incoming 
field to give rise to a signal at one port alone. 



In a previous paper [18], we developed a gen- 
eral theory for analysing phased arrays that al- 
lows overlapping beams to be modelled. The 
theory identifies the natural modes of a phased 
array, namely those orthogonal field patterns to 
which a phased array can couple, which will not 
in general be equivalent to the beam patterns. 
We showed how, once the modes have been iden- 
tified, sources can be recovered from measured 
correlations, and noise can be incorporated into 
the model. The purpose of this paper is to 
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develop the theory, illustrated by simulations, 
for a particular class of phased array: a regu- 
larly spaced array of identical, non-interacting, 
paraxial horns, equipped with a Butler matrix 
|12j to form the synthesized beams. The theory 
we developed previously is much more general, 
allowing for interaction between the horns, and 
for wide angle primary beams, as well as plac- 
ing no restriction, other than linearity, on the 
type of beam combiner. In particular the theory 
permits the modelling of phased arrays of planar 
antennas, operating at short wavelengths. Addi- 
tionally, extensive work has been done on reduc- 
ing sidelobes by the use of non-uniform array 
configurations [HI [10]. Nevertheless, we restrict 
ourselves in this paper to the simplest type of 
phased array for two reasons. Firstly, the effects 
of overlapping beams on image recovery have not 
previously been fully appreciated, and we wish 
to demonstrate the effects in the simplest possi- 
ble system. Secondly, although sidelobes can be 
reduced by the use of non-uniform arrays, there 
is still a well known tradeoff between beam width 
and sidelobe levels, so our results have implica- 
tion for more sophisticated array geometries. 

2 Theory 

Consider an array of M non-interacting horns. 
Let Q, be the unit vector pointing in the direc- 
tion of arrival (DoA) for each point on the sky, 
measured from a reference point in the aperture 
plane. Because SI is constrained to be of unit 
length, it has only 2 degrees of freedom, and the 
X and y components and are the angu- 
lar coordinates in radians of the associated point 
on the sky (throughout the paper we use super- 
scripts to denote the Cartesian components of a 
vector, and subscripts to index horns and ports). 



Let drj be an infinitesimal patch of solid angle 
with DoA vector $7, and let ^ be a patch of solid 
angle large enough to contain the primary beam 
(the beam of the horn). The signal at the mth 
horn is then given by 

y^= [ ^*^{^)E{^)dn , (1) 

where (/)m(ri) is the beam of the mth horn and 
E(ft) is the magnitude of the component of the 
electric field arriving from direction ft. In this 
paper we work in terms of scalar fields, for sim- 
plicity, but in our previous paper [18] we devel- 
oped the theory in terms of vector fields, thus 
allowing polarization to be incorporated. 

A phased array is formed by taking linear com- 
binations of the signals from the horns. Assum- 
ing that there is no scattering from the beam 
forming network into the horns, or that any such 
scattering is linear, the signal at the pth port is 

M 

Zp = ^ ^ Spmy-m 5 p = 1 ^ P , (2) 

m=l 

where Spm is the scattering matrix of the beam 
forming network (which need not be square) , and 
P is the number of ports. We can thus define a 
synthesized beam '0p(ri) associated with the pth 
port by 

M 

V'p(f2) = , (3) 

m=l 

for then we can write the output of the pth port 
as 

zp= [ rpmE{n)dn (4) 

J A 

For partially coherent sources, we can write a 
similar result in terms of coherence functions: 

(^p4)= / / V';(f^i)^p'(n2)r(ni,n2)dJ^idi72 , 

J A J A 

(5) 
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where 



(6) 



is the mutual coherence function of the field 
on the sky. For most astronomical sources, 
r(rii,ri2) win be diagonal; however if the ob- 
ject being imaged is a maser, it will be spatially 
coherent, while if a phased array is placed in the 
focal plane of a telescope, even an incoherent 
source can give rise to coherent fields. 

2.1 Image recovery 

Suppose we have measured the outputs Zp of a 
phased array and want to reconstruct the field on 
the sky. In our previous paper |18J we introduced 
the idea of using dual beams for this purpose: we 
now describe this method. 

We showed that the beams ipp{fl) can be de- 
composed thus: 



N 



(7) 



where cj„ is a real number greater than zero, N 
is a finite integer, less than or equal to both the 
number of horns M and the number of ports 
P, and Upn and V^(ri) satisfy the orthogonality 
relations 



A 



pn' 



(8) 
(9) 



We have called Vn'{^) the "eigenfields" in pre- 
vious publications: they are an orthogonal and 
minimal set of functions the phased array can 
couple to, and can thus be considered to be its 
natural modes. In the parlance of the theory of 



continuous functions, ([7| is known as a Hilbert- 
Schmidt decomposition. When the beams are 
discretized for numerical work, ([7| becomes a 
singular value decomposition (SVD). The SVD 
has been used extensively in many applications, 
and there are numerous algorithms for its com- 
putation. Its two primary uses are in data com- 
pression and the calculation of generalized in- 
verses; in this paper we use it in the latter ca- 
pacity. We showed that the dual beams ipp^ft), 
analagous to the generalized inverse of a matrix 
[T3l [Ml IT5] may be defined as 



N 



(10) 



The incident field may then be estimated by 
p 



E{n) = ZptPpiQ) . 



(11) 



(11) gives the best estimate of the incoming field 



in the sense that the error E{fl) — E{fl) is math- 
ematically orthogonal to all of the eigenfields, 
the synthesized beams, and the dual beams. In 
the case that the beams span the incoming field 
then E(fl) = E(fl); if the beams span all possi- 
ble incoming fields then we say that the beams 
constitute a frame [5] . 

In the case of the partially coherent field (|5| , 
the coherence of the incoming field can be esti- 
mated by 



p 

E 

p,p' 



)^p{n)%^{n) . (12) 



This is also a "best estimate" in the sense that 
the error T{fti,ft2) — T{fti,ft2) is orthogonal 
to any dyadic formed from the eigenfields. It is 
not however constrained to be incoherent. We 
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have addressed this problem in [18j, introduc- 
ing the idea of "interferometric frames" , and we 
will develop the idea further in future publica- 
tions; however we note that setting fti = fl2 



in ( 12 ) still provides a good estimate of the in- 



tensity distribution of the source, and indeed is 
numerically more efficient than the use of inter- 
ferometric frames. 

2.2 Array of identical paraxial horns 
with a Butler matrix 

If the horns are identical, and their beams are 
sufficiently narrow to be described by paraxial 
optics, their beams can be trivially related to the 
beam of a single horn at an arbitrary reference 
position by 



where 6^ is the displacement vector of the mth 
horn in wavelengths from the reference position. 
The normal choice of beam forming network is a 
Butler matrix, which can be represented by the 
scattering matrix 



pm 



(14) 



and gives rise to the synthesized beams 

M 

V'p(n) = e-^^™-("''~"V(fi) , (15) 

m=l 

where cxp is a DoA unit vector chosen to 
move the synthesized beams around on the sky. 
Henceforth in this paper we consider only arrays 
for which = and hence 17^ and Sp can be 
ignored. In practice a Butler matrix can be con- 
structed from a lossless combination of 90° hy- 
brid couplers and phase shifters \12\ ■ The reason 



for using a Butler matrix is that the signal from 
the pth port is 



M 



A 



e'b"-<'^p~^)^*{^)E{fl)dn , (16) 



m=l 



and, for an array with a sufficient number of 
horns, 

Zp^<f>*{ap)E{ap). (17) 

Thus, in the limit of perfect coverage, the sig- 
nal from each port corresponds to a point on the 
sky, weighted by the conjugate of the primary 
beam. However for a real array, the synthesized 
beams have finite extent, so even with a simple 
Butler matrix it is straightforward to generate 
overlapping, non-orthogonal beams; indeed it is 
necessary to choose the spacing of the horns care- 
fully if near-orthogonal beams are desired. We 



(13) note that (16) is invariant under the simultane- 



ous transformation 

n 

a 
b 



cn 

Co. 

1 



c 



b, 



(18) 



for some non-zero constant C. Thus the same 
behaviour is shown if we scale the baseline vec- 
tors in inverse proportion to the angular scales. 

We now consider the detailed behaviour of the 
synthesized beams when the coverage is not per- 
fect, in the sense that the Fourier components of 
the incident waves are not fully sampled. First 
we consider the case of a densely packed square 
array, with dimensions x wavelengths. In 



this case the summation in (15) can be approx- 



imated by an integral, resulting in the approxi- 
mation 

i;p{tt) ~ sinc(S^(f7^-S^)/2)sinc(S^(O^-ap/2)0(n) 

(19) 
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Thus the spatial frequency of the sidelobes and 
the width of the central lobe of the beam are 
determined by the size of the array. 

Now we consider an array that is less densely 
packed. The beam patterns will clearly be more 
complex, but the width of the central lobe and 
the typical spacing of the sidelobes are expected 
still to be determined largely by the overall array 
size. Another effect of the sampling however is 
to introduce grating lobes. If the position vector 
of the pth horn can be written in the form 



bp = nlAb^ei + n^Aft^es 



(20) 



where and are integers and {61,62} is a 



set of lattice vectors, then from (15) and (20) we 
see that 



<A(n + mi|^ei+m2^e2) 



We have discussed the properties of the syn- 
thesized beams; now we discuss practical issues 
relating to their use. In the past, the approach 
has been to make the synthesized beams near- 
orthogonal, but the effect of this is to leave large 
gaps on the sky that are not covered by any 
beam, and additionally to increase the level of 
the sidelobes. As we shall demonstrate, this is an 
unnecessary restriction on the design of phased 
arrays, and indeed it is desirable that overlap- 
ping beams should be used. 

3 Numerical results and discus- 
sion 

In this section we simulate phased arrays con- 
sisting of a square array of identical non- 
interacting horns, with Gaussian beam patterns, 
and equipped with Butler matrices as described 



in Section 2.2 the beams being given by (15). 



reciprocal basis satisfying • ej = 5ij. The pe- 
riodicity depends thus on the minimum spacing 
between horns, not on the overall size of the ar- 
ray, and furthermore does not require the grid 
of horns to be square or indeed for all positions 
on the grid to be filled. If the position vectors of 



for all integers and m^, where {ei, 63} is the ^e use the rescaUng given in pi). The phase 



the horns cannot be written in the form (20), the 



strict periodicity will not be seen; nevertheless, 
significant sidelobes are still likely. 

To summarize, we expect to see the following 
features in the synthesized beams of a regular 
square array of size B wavelengths and spacing 
AB wavelengths: 

1 . A principal lobe of width ~ ^ 

2. Sidelobes with typical spacing ~ ^ 

3. Grating lobes with spacing ~ 



slopes Sp are chosen such that the centers of the 
beams form a square grid of dimension C radi- 
ans, while the grid used for the numerical inte- 
grations has dimension 2C radians, and 200x200 
grid points (see Figure [T]). The horns have Gaus- 
sian beams with width 0.5C radians. The left 
hand panel of Figure [2] shows the beam of the 
horns. The dual beams were obtained from ([7]) 
and ([10|, using the LAPACK [J "divide and 
conquer" routine to perform the SVD. In cal- 
culating the duals, singular values below one 
part in 10^ of the largest singular value were 
discarded, corresponding to a dynamic range of 
lOOdB. In addition to finding the beams, duals 
and singular value decomposition, we tested the 
ability of the array to recover both a fully coher- 
ent and a fully incoherent source, using equation 



(11) and (12) respectively. The test source used 
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is shown in the right hand panel of Figure |2j 

Figures [3] and |4] show the beams and dual 
beams of a phased array with a square grid of 
5x5 horns, of overall size 10/C wavelengths, 
and with 5x5 synthesized beams such that 
P = M. The key features in the synthesized 
beams are the central maximum, which moves 
across the sky as Sp is varied, and the large 
sidelobes. As would be expected from the form 
of the beams, they are far from orthogonal, and 
the dual beams take a very different form to the 
beams. Figures |5] and [6] show the corresponding 
eigenfields (which are the natural modes) and 
singular values; the singular values are plotted 
in decibels, normalized to the largest singular 
value. There are fewer eigenfields than synthe- 
sized beams because some of the singular val- 
ues are below the dynamic range of lOOdB be- 
low the largest singular value. Note that the 
lower order eigenfields bear some resemblance to 
prolate spheroidal wavefunctions, but the higher 
order ones exhibit unusual and complex geome- 
tries. Note also that the integral operator de- 
scribing the phased array is poorly conditioned 
but non-singular, because the number of horns 
is equal to the number of ports. In Figure [7] 
we show reconstructed fully and mutually co- 
herent and incoherent fields obtained from the 
array. They have similar features, except that 
the artefacts are more prominent in the coher- 
ent case. We observe the physical significance 
of the dual beams: their geometry characterizes 
the artefacts obtained when a reconstruction of 
the field is performed using the SVD. 

The synthesized beams depend on both the 
geometry of the array and the form of the pri- 
mary beams. To separate these effects we show 
the beams, dual beams, eigenfields and spectrum 
of the same phased array as above but with a 
uniform primary beam. Clearly this is not phys- 



ically realistic, but it serves to demonstrate the 
role of the primary beam. The results are shown 
in Figures [8) [9l 10 and 11 When we compare 



these results with those obtained with a Gaus- 
sian primary beam, we observe that the Gaus- 
sian beam has the effect of truncating the syn- 
thesized beams at the edges, as expected. An 
effect of this is to shift the maxima in the sythe- 
sized beams towards the center, i.e. the maxi- 
mum of the pth synthesized beams is not in fact 
at the vector Sp. The dual beams take a much 
simpler form without the primary beam; indeed 
they differ from each other only in the phases, 
as the amplitudes are all the same. The spec- 
trum shows that the matrix is better conditioned 
when the primary beam is uniform; this is not 
surprising as information is expected to be lost 
when the beams are tapered. The lower order 
singular vectors are of similar form whether the 
primary beam is uniform or not, and are merely 
tapered at the edges by the primary beam; the 
higher order singular vectors are quite different. 
This would be expected because, for a tapered 
beam, the higher order singular vectors are con- 
centrated away from the center. 

3.1 Variable number of horns, fixed 
array size 

Now we consider the effect of changing the num- 
ber of horns, while keeping the size of the ar- 
ray and number of ports fixed. We revert to 



the Gaussian primary beam. Figures 12 and 13 
show the synthesized beams and spectrum for 
an array with 16 horns in a 4 x 4 grid (for which 
P > M); we see that the spectrum is truncated 
to 16 singular values. This is a different effect 
to the poor conditioning resulting from the over- 
lapping beams seen previously, for here the spec- 
trum falls strictly to zero. Thus the number of 
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horns places an absolute restriction on the num- 
ber of degrees of freedom of the phased array, 
while the beam overlap imposes a restriction de- 
pending on the dynamic range, which in turn 
depends on a number of factors including the 
noise levels. The effect of the reduced number of 
horns is manifest in the increased magnitude and 
higher spatial frequency of the sidelobes of the 
synthesized beams. The higher spatial frequency 
might be expected from the fact that the gaps 
between horns are larger, but we observe that 
this is only evident when the number of horns is 
fewer than the number of ports. 

We also increased the number of horns to 
10 X 10 (such that P < M), and found that 
the results were identical to those with 5x5 
ports (data not shown). This can be explained 
by the fact that the maximum number of degrees 
of freedom is dictated by the smaller of the num- 
ber of ports and the number of horns; there is 
therefore no advantage in increasing the number 
of horns above the number of ports. 

3.2 Variable array size, fixed number 
of horns 

Now we consider the effect of changing the size 
of the array, while maintaining the same num- 
ber of horns. We revert to the original number 
of horns, in a 5 x 5 array, but with the array 



dimensions halved. Figures 14 , 15 and 16 show 



the beams, duals and spectrum of such an ar- 
ray. The size of the synthesized beams increases, 
while the spatial frequency of sidelobes reduces 
so much that they are outside the primary beam, 
as would be expected. The dual beams take on a 
more complicated and less symmetric form. The 
spectrum falls more rapidly, as expected from 
the strongly overlapping beams. Remarkably, it 
is still possible to recover images with such an 



array, as shown in Figure 17 Again, the recon- 
structed fields contain artefacts resembling the 
dual beams. In practise the quality of image re- 
covery is limited by the dynamic range of the 



system, as we shall show in Section 3.3 



Now we increase the size of the array, again 
keeping the number of horns constant. Figures 



18 and 19 shows the synthesized beams and du- 



als for an array of twice the original dimensions, 
while Figure 20 shows the spectrum. We see 
that as the array size is increased the beams 
get smaller, and the frequency of the sidelobes 
increases. The beams also become more or- 
thogonal, such that the duals closely resemble 
the beams, and the conditioning also improves. 
However, the reconstructed images in Figure 21 
are of poor quality, and appear to be inferior 
to those obtained with the strongly overlapping 



beams in Figure 17 The reason is clear: in or- 
der to make the central lobes of the synthesized 
beams narrow enough that the beams are orthog- 
onal, it is necessary to increase the frequency of 
the sidelobes, thus introducing artefacts into the 
reconstructed field. Moreover, the higher the fre- 
quency of the the sidelobes, the greater their in- 
tensity relative to the main lobe, and hence the 
greater the difficulty in deconvolving the image 
will be. 



Figures 22, 23 and 24 show the beams, duals 
and spectrum of an array of 4 times the orig- 
inal dimensions, and the same effects are seen 



as in Figures 18 -20 only more exaggerated. In 
addition however we see the periodicity in the 
synthesized beams discussed in Section |2.2[ and 
furthermore the principle maxima in the beams 
are located in the middle, not at the "correct", 
positions, as a result of truncation by the pri- 
mary beam. This leads to aliasing, and hence to 
the completely erroneous reconstructed fields in 
Figure |25l 
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3.3 Dynamic range 

The threshold below which singular values must 
be discarded, which dicates the value of A'' in ([T]) , 
is determined by the dynamic range of the sys- 
tem. In the preceding simulations we have used 
an implied dynamic range of lOOdB; we now con- 
sider the effect of reducing this to 50dB. In Fig- 
ure [26] we show the reconstructed fields obtained 
with this reduced dynamic range for the array of 



dimensions 5/C wavelengths used in Figure 17 



We see that the reconstructed coherent field is of 
inferior quality with the reduced dynamic range. 
At first sight, the reconstructed incoherent field 
appears not to be so bad, but on closer inspec- 
tion we see that the maxima have been shifted 
towards the center of the image, and away from 
their true position. This is because the discarded 
eigenfields are concentrated away from the cen- 
ter and thus contain the necessary information 
to correctly determine the position of the source. 
The dynamic range of the instrument is thus cru- 
cially important in determining the quality of 
the reconstructed image. The question we now 
ask is whether, in the case of limited dynamic 
range, we can improve the situation by synthe- 
sizing more beams. 

3.4 Increasing number of synthesized 
beams 

We now examine whether increasing the number 
of synthesized beams for a given configuration of 
horns makes any difference to the ability of the 
array to recover a field. We performed simula- 
tions for an array of dimension 5/C, like that 



used in Figures 14 -17 but with 100 synthesized 
beams in a 10 x 10 grid covering the same patch 
of sky. The results were unchanged to Figures 
[l4|jl7| (data not shown). Even when the dynamic 



range is reduced to 50dB, there is no difference 
obtained by using more beams. 

We also consider whether increasing the num- 
ber of beams has an effect for the array of dimen- 



sion 20/C, as used in Figure 21 , In this case, the 



rank is equal to the number of horns; the ques- 
tion is whether the sidelobes can be removed by 
synthesizing more beams. Increasing the number 
of ports to 100 once again made no difference; 
our results strongly suggest that there is no ad- 
vantage, at least in terms of modal throughput, 
in making the number of ports different to the 
number of horns. 



4 Conclusion 

We have performed simulations of imaging 
phased arrays, using our previously developed 
theory to calculate eigenfields and dual beams, 
and to reconstruct both coherent and incoherent 
test sources. The beam patterns show all the 
features one would expect for beams produced 
by a Butler matrix; the dual beams show com- 
plicated behaviour that results in artefacts in the 
reconstructed sources. 

We have demonstrated that mathematically 
orthogonal beams do not give the best image 
quality, due to the inevitably increased sidelobe 
frequency. On the contrary some degree of over- 
lap is desirable; however the overlap should ide- 
ally not be so much that the range of singular 
values exceeds the dynamic range of the system, 
as this degrades the image. We have looked for 
and failed to find any advantage in making the 
number of horns different to the number of ports. 

The significance for the design of imaging 
phased arrays seems to be that, as a general rule, 
the number of ports should equal the number of 
horns, and that the beams should be made to 
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overlap as much as the dynamic range of the in- 
strument permits. It follows that improving the 
dynamic range of a phased array, by reducing 
noise for instance, will improve the quality of 
the images not only by reducing the noise in the 
image, but also by reducing the artefacts. 
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Figure 3: Linear grayscale plots showing the 
magnitude, as a function of Jl^'/C and ^I'^/C, 
of the synthesized beams of an array with 5x5 
horns, total size 10 x 10/C wavelengths 



Figure 1: Positions of the sythesized beams on 
the sky 
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Figure 2: Linear grayscale plots showing the 
magnitude, as a function of fi^/C and of 
the primary beam (left) and test source (right) 
used in the subsequent simulations 
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Figure 4: Linear grayscale plots showing the 
magnitude, as a function of fl^/C and Q^/C, 
of the dual beams of an array with 5x5 horns, 
total size 10 x 10/C wavelengths 
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Figure 7: Linear grayscale plots showing the 
magnitude of a reconstructed coherent source 
(left), and the intensity of a reconstructed in- 
coherent source (right), as a function of 
and ri^/C, for an array with 5x5 horns, total 
size 10 X 10/C wavelengths 



Figure 5: Linear grayscale plots showing the 
magnitude, as a function of Q^/C and W /C, 
of the eigenfields of an array with 5x5 horns, 
total size 10 x 10/C wavelengths 
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Figure 6: Singular values in dB of an array 
with5 X 5 horns, total size 10 x 10/C wavelengths 
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Figure 8: Linear grayscale plots showing the 

magnitude, as a function of fi^/C and O^/C, 
of the synthesized beams of an array with 5x5 
horns, total size 10 x 10/C wavelengths, and uni- 
form primary beams 
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Figure 9: Linear grayscale plots showing the 
magnitude, as a function of ^l^/C and ^I'^/C, 
of the dual beams of an array with 5x5 horns, 
total size 10 x 10/C wavelengths, and uniform 
primary beams 
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Figure 10: Linear grayscale plots showing the 
magnitude, as a function of VL^ /C and W> jC ^ of 
the eigenfields of an array with 5x5 horns, total 
size 10 X 10/C wavelengths, and uniform primary 
beams 



5 to 13 2D 25 



Figure 11: Singular values in dB of an array 
with5 X 5 horns, total size 10 x 10/C wavelengths, 
and uniform primary beams 
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Figure 12: Linear grayscale plots showing the 
magnitude, as a function oifl^/C and W /C, of 
the synthesized beams of an array with 4x4 
horns, total size 10 x 10/C wavelengths 
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Figure 14: Linear grayscale plots showing the 
magnitude, as a function of il^/C and il^'/C, of 
the synthesized beams of an array with 5x5 
horns, total size 5 x 5/C wavelengths 
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Figure 13: Singular values in dB of an array with 
4x4 horns, total size 10 x 10/C wavelengths 



Figure 15: Linear grayscale plots showing the 
magnitude, as a function of Cl^/C and ^1^/0, of 
the synthesized beams of an array with 5x5 
horns, total size 5 x 5/C wavelengths 
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Figure 16: Singular values in dB of an array with 
5x5 horns, total size 5 x 5/C wavelengths 
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Figure 17: Linear grayscale plots showing the 
magnitude of a reconstructed coherent source 
(left), and the intensity of a reconstructed in- 
coherent source (right), as a function of fl^/C 
and ri^/C, for an array with 5x5 horns, total 
size 5 X 5/C wavelengths 
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Figure 18: Linear grayscale plots showing the 
magnitude, as a function of Q^/C and W /C, of 
the synthesized beams of an array with 5x5 
horns, total size 20 x 20/C wavelengths 
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Figure 19: Linear grayscale plots showing the 
magnitude, as a function of and O^^/C, of 

the dual beams of an array with 5x5 horns, total 
size 20 X 20/C wavelengths 




Figure 20: Singular values in dB of an array with 
5x5 horns, total size 20 x 20/C wavelengths 
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Figure 21: Linear grayscale plots showing the 
magnitude of a reconstructed coherent source 
(left), and the intensity of a reconstructed in- 
coherent source (right), as a function of 
and ri^/C, for an array with 5x5 horns, total 
size 20 X 20/C wavelengths 
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Figure 22: Linear grayscale plots showing the 
magnitude, as a function of i7^/C and 0,^/0, of 
the synthesized beams of an array with 5x5 
horns, total size 40 x 40/C wavelengths 




Figure 23: Linear grayscale plots showing the 
magnitude, as a function of /C and W^/C, of 
the dual beams of an array with 5x5 horns, total 
size 40 X 40/C wavelengths: duals 



Figure 24: Singular values in dB of an array with 
5x5 horns, total size 40 x 40/C wavelengths 



Figure 25: Linear grayscale plots showing the 
magnitude of a reconstructed coherent source 
(left), and the intensity of a reconstructed in- 
coherent source (right), as a function of Q^/C 
and 0^/C, for an array with 5x5 horns, total 
size 40 X 40/C wavelengths 
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Figure 26: Linear grayscale plots showing the 
magnitude of a reconstructed coherent source 
(left), and the intensity of a reconstructed inco- 
herent source (right), as a function oiCl^/C and 
ri^/C, for an array with 5x5 horns, total size 
5 X 5/C wavelengths, and dynamic range 50dB 
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